Fingerprints of spin-fermion pairing in cuprates 
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We demonstrate that the feedback effect from bosonic excitations on fermions, which in the past 
allowed one to verify the phononic mechanism of a conventional, s— wave superconductivity, may 
also allow one to experimentally detect the "fingerprints" of the pairing mechanism in cuprates. We 
argue that for spin-mediated d— wave superconductivity, the fermionic spectral function, the density 
of states, the tunneling conductance through an insulating junction, and the optical conductivity are 
affected by the interaction with collective spin excitations, which below T c are propagating, magnon- 
like quasiparticles with gap A 3 . We show that the interaction with a propagating spin excitation 
gives rise to singularities at frequencies A + A 3 for the spectral function and the density of states, 
and at 2A + A s for tunneling and optical conductivities, where A is the maximum value of the 
d— wave gap. We further argue that recent optical measurements also allow one to detect subleading 
singularities at 4A and 2 A + 2A S . We consider the experimental detection of these singularities as 
a strong evidence in favor of the magnetic scenario for superconductivity in cuprates. 



PACS numbers:71.10.Ca,74.20.Fg,74.25.-q 
I. INTRODUCTION 

One of the very few accepted facts for high-T c materi- 
als is that they are d-wave superconductors Jl| ||. This 
salient universal property of all cuprates entails strong 
constraints on the microscopic mechanism of supercon- 
ductivity. However, it does not uniquely determine it, 
leading to a quest for experiments which can identify 
"fingerprints" of a specific microscopic theory of rf-wave 
superconductivity, a strategy somewhat similar to the 
one used in conventional superconductors (see e.g., [Q). 
There, the identification of characteristic phonon fre- 
quencies in the tunneling density of states (DOS) below 
T c was considered as a decisive evidence for the electron- 
phonon mechanism for superconductivity [^j. 

In this paper, we assume a'priori that the pairing in 
cuprates is mediated by the exchange of collective spin 
excitations. It has been demonstrated both numerically 
and analytically that this exchange gives rise to a d-wave 
superconductivity Q. We discuss to which extent the 
"fingerprints" of the spin-mediated pairing can be ex- 
tracted from the experiments on hole-doped high T c ma- 
terials. We argue that due to strong spin-fermion cou- 
pling, there is a very strong feedback from spin excita- 
tions on fermions, specific to d— wave superconductors 
with magnetic pairing interaction. The origin of this 
feedback is the emergence of a propagating collective 
spin bosonic mode below T c . We show that this mode 
is present for any coupling strength, and its gap A s is 
smaller than the minimum energy ~ 2A which is neces- 
sary to break a Cooper pair. In the vicinity to the anti- 
ferromagnetic phase, A s cx £~ 1 where £ is the magnetic 
correlation length. We show that the spin propagating 
mode changes the onset frequency for single particle scat- 
tering, and gives rise to the "peak-dip-hump" features in 
angular resolved photoemission (ARPES) experiments, 



the "dip-peak" features in tunneling experiments, and to 
the singularities and fine structures in the optical conduc- 
tivity. We demonstrate that (i) these features have been 
observed (ii) ARPES R,B ,|l2|,p| , tunnel- 

ing Jl^ , |l5| , and conductivity data jl6|L8| are consistent 
with each other, and (iii) the value of A s extracted from 
these various experiments agrees well with the resonance 
frequency measured directly in neutron scattering exper- 
iments 



A. The physical origin of the effect 

The physical effect which accounts for dips and humps 
in the density of states and spectral function of cuprates 
by itself is not new and is known for conventional s— wave 
superconductors as a Holstein effect [p|,|22|J23|] 





FIG. 1. a) The exchange diagram for boson mediated in- 
teraction. The solid line stands for a propagating fermion. 
The wiggled line is a phonon propagator in case of electron- 
phonon interaction, and a magnon line in case of spin- fluctu- 
ation mediated interaction, b) The lowest order diagram for 
the fermionic self energy due to a direct four fermion interac- 
tion, also represented by a wiggly line 

Consider a clean s— wave superconductor, and suppose 
that the residual interaction between fermions occurs via 
the exchange of an Einstein phonon. Assume for simplic- 
ity that the fully renormalized electron phonon coupling 
is some constant g ep , and that the phonon propagator 
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D(q,w) is independent on the momentum q and has a 
single pole at a phonon frequency, A p (a Holstein model) 
p2| pif . The phonon exchange gives rise to a fermionic 
self-energy (see Fig |l]a) 



£(w m ) 



lLO r . 



(1) 



which is a convolution of -D(w) = 1/ (A 2 — (u> + iS) ) with 
the full fermionic propagator Gk{to). (For further conve- 
nience, we absorbed a bare iu> m term into the definition of 
S(w m )). In a superconductor, the fermionic propagator 
is given by 



G k (w) 



S(w) +£ k 



E 2 (w)-$ 2 (w)-£ 2 



(2) 



where $(w) is the anomalous vertex function, and £k 
is the band dispersion of the fermions. The supercon- 
ducting gap, introduced in the BCS theory, is the solu- 
tion of S(A) = $(A). [Alternatively to and $(w), 
one can introduce the complex effective mass function, 
Z(lo) — E(w)/w, and the complex effective gap function 
A (a;) = $(J)/Z{u)) In what follows we will use 
£'(w), &(<jj), etc. to denote real parts and £"(w), $"(w), 
etc. for the imaginary parts of the functions we study. 

For T = one can rigorously prove that both £"(w) 
and $"(w) vanish for W < A. This implies that the 
fermionic spectral function A^(uj) = \G'^(ui)\ /it for par- 
ticles at the Fermi surface (k = kp) has a 6— functional 
peak at u) = A, i.e. A is a sharp gap in the excitation 
spectrum at zero temperature. Also, the fermionic den- 
sity of states in a superconductor 



N(oj) = Im 



£(w) 



($ 2 (w) - E 2 (o))) 1 / 2 



(3) 



vanishes for a; < A and has a square-root singularity 
N(ll>) oc (w — A) -1 / 2 for frequencies above the gap, to > 
A. 

The onset of the imaginary part of the self-energy due 
to inelastic single particle scattering can be easily ob- 
tained by applying the spectral representation to Eq. (|l|) 
and re-expressing the momentum integration in terms of 
an integration over £k- At T = we then obtain 



£"(w > 0) cx / duj'N{u)')D"(u) - u>') 



(4) 



Since for positive frequencies, D"(lo) — (ttDo/2A p )6(uj — 
A p ), the frequency integration is elementary and yields 



E"(lo > 0) cx N(u-Ap). 



(5) 



We see that the single particle scattering rate is di- 
rectly proportional to the density of states shifted by 
the phonon frequency. Clearly, the imaginary part of 



the fermionic self-energy emerges only when u> exceeds a 
threshold at 



w 



A + A„, 



(6) 



i.e., the sum of the superconducting gap and the phonon 
frequency. Right above this threshold, £"(w) oc (ui — 



u ) 



~ x / 2 . By Kramers-Kronig relation, this nonanalytic- 
ity causes an analogous square root divergence of S'(w) 
at lj < wo- Combining the two results, we find that near 
the threshold, E(w) = A + C / y/uja — u where A and C 
are real numbers. By the same reasons, the anomalous 
vertex $(w) also possesses a square-root singularity at 
loq. Near lu — loq, = B + C/y/uo — uj with real B. 

Since u>o > A, we have A > B. 

The singularity in the fermionic self-energy gives rise 
to an extra dip-hump structure of the fermionic spectral 
function at k = k^?. Below loq, the spectral function is 
zero except for w = A, where it has a 5— functional peak. 
Above w , A(lu) cxIm(E(w)/(E 2 (w)-$ 2 (w)) emerges as 
A{uj) cx (oj — wo) 1 / 2 . At larger frequencies, A(uj) passes 
through a maximum, and eventually vanishes. Adding a 
small damping due to either impurities or finite temper- 
atures, one obtains the spectral function with a peak at 
lu = A, a dip at w loq, and a hump at a somewhat 
larger frequency. This behavior is schematically shown 
in Fig. §. 




FIG. 2. The schematic form of the quasiparticle spectral 
function in an s— wave superconductor. Solid line - T — 0, 
dashed line - a finite T. u>o — A + A p 

The singularities in S(w) and <&(w) affect other ob- 
servables such as fermionic DOS, optical conductivity, 
Raman response, and the SIS tunneling dynamical con- 
ductance p2|p7| . 

For a more complex phonon propagator, which de- 
pends on both frequency and momentum, actual singu- 
larities in the fermionic self-energy and other observables 
are weaker and may only show up in the derivatives over 
frequency [ p6[ . Still, however, the opening of the new re- 
laxational channel at wo gives rise to singularities in the 
electronic properties of an s— wave superconductor. 
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B. The similarities and discrepancies between d- and 
s— wave superconductors 

For magnetically mediated d— wave superconductivity, 
the role of phonons is played by spin fluctuations. As 
we said, these excitations are propagating, magnon-like 
modes below T c (more accurately, below the onset tem- 
perature for the pseudogap), with the gap A s . This spin 
gap obviously plays the same role as A p for phonons, 
and hence we expect that the spectral function should 
display a peak-dip-hump structure as well. We will 
also demonstrate below that for the observables such 
as the DOS, Raman intensity and the optical conduc- 
tivity, which measure the response averaged over the 
Fermi surface, the angular dependence of the d— wave 
gap A(6) cx cos (29) softens the singularities, but does 
not wash them out over a finite frequency range. Indeed, 
we find that the positions of the singularities are not de- 
termined by some averaged gap amplitude but by the 
maximum value of the d— wave gap, A(0) = A. 

Despite similarities, the feedback effects for phonon- 
mediated s— wave superconductors, and magnetically 
mediated d— wave superconductors arc not equivalent as 
we now demonstrate. The point is that for s— wave super- 
conductors, the exchange process shown in Figjlja is not 
the only possible source for the fermionic decay: there ex- 
ists another process, shown in Fig.^Jb, in which a fermion 
decays into three other fermions. This process is due to a 
residual four-fermion interaction |^,^7|. One can easily 
make sure that this second process also gives rise to the 
fermionic decay when the external lu exceeds a minimum 
energy of 3A, necessary to pull all three intermediate 
particles out of the condensate of Cooper pairs. At the 
threshold, the fermionic spectral function is non-analytic, 
much like at A + A p . This implies that in s-wave super- 
conductors, there are two physically distinct singularities, 
at A + A p and at 3 A, which come from different processes 
and therefore are independent of each other. Which of 
the two threshold frequencies is larger depends on the 
strength of the coupling and on the shape of the phonon 
density of states. At weak coupling, A p is exponentially 
larger than A, hence 3 A threshold comes first. At strong 
coupling, A s and A are comparable, but calculations 
within Eliashberg formalism show that for real materials 
( e.g. for lead or niobium) still, 3A < A + A p . g]. This 
result is fully consistent with the photoemission data for 
these materials [[L7). 

For magnetically mediated d-wave superconductors the 
situation is different. In the one-band model for cuprates, 
which we adopt, the underlying interaction is solely a 
Hubbard-type four-fermion interaction. The introduc- 
tion of a spin fluctuation as an extra degree of freedom is 
just a way to account for the fact that there exists a par- 
ticular interaction channel, where the effective interac- 
tion between fermions is the strongest due to a closeness 



to a magnetic instability. This implies that the propaga- 
tor of spin fluctuations by itself is made out of particle- 
hole bubbles like those in Fig.|l|b. Then, to the lowest 
order in the interaction, the fermionic self-energy is given 
by the diagram in Fig.^Jb. Higher-order terms convert a 
particle-hole bubble in Fig.|l]b. into a wiggly line, and 
transform this diagram into the one in Fig.|l|a. Clearly 
then, a simultaneous inclusion of both diagrams would 
be a double counting, i.e., there is only a single process 
which gives rise to the threshold in the fermionic self- 
energy. 

Leaving a detailed justification of the spin-fermion 
model to the next section, we merely notice here that 
the very fact that the diagram in Fig [l|b is a part of that 
in Fig.pi implies that the development of a singularity 
in the spectral function at a frequency different from 3A 
cannot be due to effects outside the spin-fermion model. 
Indeed, we will show that the model itself generates two 
singularities, at 3 A and at A + A s < 3 A. The fact 
that this is an internal effect, however, implies that A s 
depends on A. The experimental verification of this de- 
pendence can then be considered as a "fingerprint" of the 
spin-fluctuation mechanism. Furthermore, as the singu- 
larities at 3 A and A + A s are due to the same interaction, 
their relative intensity is another gauge of the magnetic 
mechanism for the pairing. We will argue below that 
some experiments on cuprates, particularly the measure- 
ments of optical conductivity ]lS| ], allow one to detect 
both singularities, and that their calculated relative in- 
tensity is consistent with the data. 

II. SPIN-FERMION MODEL 

The point of departure for our analysis is the effec- 
tive low-energy theory for Hubbard-type lattice fermion 
models. As mentioned above, we assume a'priori that 
integrating out states with high fermionic energies in a 
renormalization group sense, one obtains low-energy col- 
lective bosonic modes only in the spin channel. In this 
situation, the low-energy theory should include fermions, 
their collective bosonic spin excitations, and the interac- 
tion between fermions and spins. This model is called a 
spin-fermion model p9| , and is described by 

n = J2 v F-( k - k F)cl a c Ka + J2 xo 1 (q)s q S- q + 

k.a q 

q,k,a,/3 

Here is the fermionic creation operator for an elec- 
tron with crystal momentum k and spin a, Oi are the 
Pauli matrices, and g is the coupling constant which mea- 
sures the strength of the interaction between fermions 
and their collective bosonic spin degrees of freedom, 
characterized by the spin-1 boson field, S q . The latter 
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are characterized by a bare spin susceptibility xo(q) = 
Xo£ 2 /(l + (q— Q) 2 £ 2 )i where £ is the magnetic correlation 
length. 

The relevant topological variable in the theory is the 
shape of the Fermi surface. We assume that the Fermi- 
surface is hole- like, i.e., it is centered at (it, it) rather than 
at (0,0). This Fermi surface is consistent with the pho- 
toemission measurements for Bi2212, at least at and be- 
low optimal doping. Luttinger theorem implies that this 
Fermi surface necessary contains hot spots - the points 
at the Fermi surface separated by the antiferromagnetic 
momentum Q. In Bi2212, these hot spots are located 
near (0, ir) and symmetry related points Jt]||. The exact 
location of hot spots is however not essential for our cal- 
culations. It is only important that hot spots do exist, 
are not located close to the nodes of the superconducting 
gap, and that the Fermi velocities at kh s and kh s + Q 
are not antiparallel to each other. We also assume that 
ujq is smaller then eo )7r , i-e., the van-Hove singularity of 
the electronic dispersion at k — (0, tt) is irrelevant for our 
analysis. 

Observe also that our bare Xo(q) does not depend on 
frequency. In general, the integration over high energy 
fermions may give rise to some frequency dependence 
of xo(Q: w )- However, as xo(q, w) comes from fermions 
with lo ~ E-p-i its frequency dependence holds in powers 
of Lo/Ep) 2 . We will see that this frequency dependence 
can be safely neglected as it is completely overshadowed 
by the iuig 2 xo/v F term which comes from low-energy 
fermions. The presence of hot spots is essential in this 
regard because a spin fluctuation with a momentum near 
Q can decay into fermions at or near hot spots. By virtue 
of energy conservation this process involves only low- 
energy fermions and therefore is fully determined within 
the model. 

This evolution of the bosonic dynamics from propagat- 
ing to relaxational is the key element which distinguishes 
between spin-mediated and phonon-mediated supercon- 
ductivities. For phonon superconductors, the bosonic 
self-energy due to spin-fcrmion interaction also contains a 
linear in to term. However, this term has an extra relative 
smallncss in u/vp cx m/M, where u is the sound velocity, 
m is the electron mass, and M is the mass of an ion. Due 
to this extra smallness, the linear in u> term in the phonon 
propagator becomes relevant only at extremely low fre- 
quencies, unessential for superconductivity. At lo ~ A, 
the bare to 2 term dominates, i.e., the phonon propaga- 
tor preserves its input form. Alternatively speaking, the 
renormalization of the phonon propagator by fermions is 
a minor effect while the same renormalization of the spin 
propagator dominates the physics below Ep. 



A. Computational technique 

The input parameters in Eq. (Q) are the coupling con- 
stant g, the spin correlation length, £, the Fermi velocity 
vp (which we assume to depend weakly on the position 
on the Fermi surface), and the overall factor xo- The lat- 
ter, however, can be absorbed into the renormalization 
of the coupling constant g — g 2 X0j an d should not be 
counted as an extra variable. Out of these parameters 
one can construct a dimensionless ratio A = "ig / Attvf^ -1 
and an overall energy scale Q = 9g/167r (the factors 3/47T 
and 9/167T are introduced for further convenience). All 
physical quantities discussed below can be expressed in 
terms of these two parameters (and the angle between 
Vk h5 and Vk h3 +Q, which does not enter the theory in any 
significant manner as long as Vk hs and Vk hB +Q are not 
antiparallel to each other). One can easily make sure 
that in two dimensions, a formal perturbation expansion 
holds in powers of A. The limit A <C 1 is perturbative 
and is probably applicable only to strongly overdoped 
cuprates . The situation in optimally doped and under- 
doped cuprates most likely corresponds to a strong cou- 
pling, A > 1. The most direct experimental indication 
for this is the absence of a sharp quasiparticle peak in 
the normal state ARPES data in materials with doping 
concentration equal to or below the optimal one @||. 

At strong coupling, a conventional perturbation theory 
does not work, but we found earlier that a controllable ex- 
pansion is still possible if one formally treats the number 
of hot spots in the Brillouin zone N = 8 as a large num- 
ber. |3(],[32],[33| The justification and a detailed descrip- 
tion of this procedure is beyond the scope of the present 
paper. We refer the reader to the original publications, 
and quote here only the result: near hot spots, one can 
obtain a set of coupled integral equations for three com- 
plex variables: the anomalous vertex $k(^) ~ ^ > k hB ( w ) 
subject to the d-wave constraint $k(w) = — $k+Q(^), 
the fermionic self-energy Ek(u>) rs Ek hB (a;), an d the spin 
polarization operator Hq(u>). The anomalous vertex and 
the fermionic self-energy are related to the normal and 
anomalous Green's functions as 



GkH 
F k (w) 



E k (w) + £ k 



(8) 
(9) 



and the polarization operator is related to the fully renor- 
malized spin susceptibility as 



x(q,w) 



i + (q - Q) 2 e - n q H/ Wsf 



(10) 



We normalized n q (w) such that in the normal state 
Hq(lu) — iuj (see below). This normalization implies 
that Wgf = Cu/AX 2 
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In Matsubara frequencies the set of the three equations 
has the form 



7TT ^ <F, 



2 V^ f + n n _ 



1/2 



7TT 



E 



n m = ttT 2^ i 



1/2 



(11) 

(12) 
(13) 



Here, $ m = $ khs (uj m ) and S m = S khs (cj m ) with 
fcrmionic Matsubara frequency oj rn = (2m + 1)ttT and 
II m = IIq(cj to ) with bosonic Matsubara frequency 
uj m = 2mirT, respectively. The first two equations are 
similar to the Eliashberg equations for conventional su- 
perconductors. The presence of the third coupled equa- 
tion for II is peculiar to the spin-fluctuation scenario, and 
reflects the fact that the spin dynamics is made out of 
fermions. As in a conventional Eliashberg formalism, the 
superconducting gap A at T = is defined as a solution 
of = after analytical continuation to the real 

axis. 

The region around a hot spot where these equations are 
valid (i.e., the "size" of a hot spot) depends on frequency 
and is given by |k - k hs | ~ ^(1 + MMf) 1/2 ©H- 
We will see that at strong coupling (A > 1), the super- 
conducting gap A ~ uj s {X 2 . Hence for frequencies com- 
parable to or larger than A, typical |k — kh s | ~ A£ 1 ~ 
kp(g / Ep). In practice, g is comparable to Ep (in the 
RPA approximation for an effective one-band Hubbard 
model for CuC>2, j«(/~2 - 3eV, while Ep is com- 
parable to a bandwidth which has the same order of 
magnitude). In this situation, the self energy and the 
anomalous vertex at the hot spots are characteristic for 
the behavior of these functions in a substantial portion 
of the Fermi surface, leading to an effective momentum 
independence of the fermionic dynamics away from the 
nodes of the gap. Of course, near zone diagonals, there is 
a different physics at low frequencies, associated with the 
fact that the superconducting gap vanishes for momenta 
along the diagonals. Below we show that the physics 
close to the nodes is universally determined by the van- 
ishing superconducting gap, and therefore insensitive to 
strong coupling effects which bear fingerprints of the pair- 
ing mechanism. For these reasons we will mostly concen- 
trate our analysis to describe peculiarities of the d-wave 
state at frequencies to > A max . 



III. THE SPIN POLARIZATION OPERATOR 

Since our goal is to find the "fingerprints" of spin ex- 
citations in the fermionic variables, we first discuss the 
general form of the spin polarization operator. We show 



that in the normal state (ignoring pseudogap effects), 
gapless fermions cause a purely diffusive spin dynamics. 
However, in a d— wave superconducting state, a gap in 
the single particle dynamics gives rise to "particle like" 
propagating magnons. 



A. Normal state spin dynamics 

In the normal state (&(u>) = 0) the polarization oper- 
ator can be computed explicitly even without the knowl- 
edge of the precise form of E(a;). Indeed, from (|l3| ) we 
immediately obtain 



II m = 7rT2j(l - sign(S, i )sign(S n+m )) 



(14) 



The result depends only on the sign of the self-energy, 
but not on its amplitude. As sign(S„) = sign(2n + 1), 
the computation is elementary, and we obtain that for 
any coupling strength 



H 



(15) 



On the real axis this translates into a purely diffusive 
and overdamped spin dynamics with H(oj) — iui. This 
result holds as long as the linearization of the fermionic 
dispersion near the Fermi surface remains valid, i.e., up to 
energies comparable to the Fermi energy. The linearity 
of n(w) and thus the fact that spin excitations in the 
normal state can propagate only in a diffusive way is a 
direct consequence of the presence of hot spots. 



B. Spin dynamics in the superconducting state 

The opening of the superconducting gap changes this 
picture. Now quasiparticles near hot spots are gapped, 
so a spin fluctuation can decay into a particle-hole pair 
only when it can pull two particles out of the conden- 
sate of Cooper pairs. This implies that the decay into 
particle hole excitations is only possible if the exter- 
nal frequency is larger than 2 A. At smaller frequen- 
cies, H"(ui) — for T — 0. This result also read- 
ily follows from Eq. (|l3|) . The Kramers-Kronig relation 
n'(w) = (2/tt) f™IL"(x)/(x 2 - u 2 ) then implies that 
due to the drop in ~n"(uj), the spin polarization operator 
in a superconductor acquires a real part, which at low 
uj is quadratic in frequency: n(w) cx uj 2 /A. An essen- 
tial point for our consideration is that H(oj = 0) = 
for any S(w) and $(w). This physically implies that the 
development of the gap does not change the magnetic 
correlation length, a result which becomes evident if one 
notices that d-wave pairing involves fermions from oppo- 
site sublattices. 

Substituting the result H(ui) oc u> 2 /A into Eq.@, we 
find that at low energies, spin excitations in a d-wave 
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superconductor are propagating, gaped magnon-like ex- 
citations: 



(16) 



The magnon gap, A s ~ (Ao^f) 1 / 2 , and the magnon ve- 



locity, 



\A/g, are entirely determined by the dy- 



namics of the fcrmionic degrees of freedom. 

Eq. ( |l6| ) is meaningful only if A s < A, i.e., ay < A. 
Otherwise the use of the quadratic form for II (a;) is not 
justified. To find how A depends on the coupling con- 
stant, one needs to carefully analyze the full set Eqn. (pd[- 
|l3| ). This analysis is rather involved ^||,^3), and is not 
directly related to the goal of this paper. We skip the 
details and just quote the result. It turns out that at 
strong coupling, A < 1, i.e. for optimally and under- 
doped cuprates, the condition A > ay is satisfied as the 
gap scales with Q and saturates at A w 0.35<Z> = 0.06g 
at A = oo, when u s f oc A~ 2 = 0. In this situation, 
the spin excitations in a superconductor are propagat- 
ing, particle- like modes with the gap A s . However, in 
distinction to phonons, these propagating magnons get 
their identity from a strong coupling feedback effect in 
the superconducting state. 

At weak coupling, the superconducting problem is of 
BCS type, and A ~ w s fexp(— A -1 ) <C ay. This result 
is intuitively obvious as ay plays the role of the Debye 
frequency in the sense that the bosonic mode which me- 
diates pairing decreases at frequencies above ay. We 
see that at weak coupling, the quadratic approximation 
for II(u;) does not lead to a pole in x(Q,o>). Still, the 
pole in x"(Q,o>) does exist also at weak coupling as we 
now demonstrate. Indeed, consider II" (w) at w « 2 A. 
One can easily make sure that at this u>, one can simul- 
taneously set both fermionic frequencies in the bubble 
to be close to A, and get a strong singularity due to 
vanishing of \/£ 2 — for both fermions. Substituting 
E 2 (w) — $ 2 (o;) cx lu — A into ( fl3| ) and using the spectral 
representation, we obtain for lu = 2A + e 



U"(uj) cx 



dx 



o (z(e-s))Va- 



(17) 



Evaluating the integral, we find that II" undergoes a fi- 
nite jump at u> = 2A. By Kramers-Kronig relation, this 
jump gives rise to a logarithmic singularity in II' (oi) at 
lu = 2A: 



n'M = - [°° dx ^fW ocAiog- 

7T J2A x - V 



2A 



2AI 



(18) 




FIG. 3. Schematic behavior of the real (dashed line) and 
imaginary (solid line) parts of the particle hole bubble in the 
superconducting state. Due to the discontinuous behavior 
of II" (w) at uj = 2A, the real part II' (u) is logarithmically 
divergent at 2A. For small lu, the real part behaves like lu 2 /A. 

The behavior of II' (w) and n"(a>) is schematically 
shown in Fig j^. The fact that II' (w) diverges logarithmi- 
cally at 2 A implies that no matter how small A/oj s f is, 
x(Q, lu) has a pole at A s < 2 A, when II" (w) is still zero. 
Simple estimates show that for weak coupling, where 
ay S> A, the singularity occurs at A s = 2A(1 — Z s ) 
where Z s oc e~ Ws/ /( 2A - ) is also the spectral weight of the 
resonance peak in this limit . 

We see therefore that the resonance in the spin suscep- 
tibility exists both at weak and at strong coupling. At 
strong coupling, the resonance frequency is A s ~ A/ A <C 
A, i.e., the resonance occurs in the frequency range where 
spin excitations behave as propagating magnon-like ex- 
citations. At weak coupling, the resonance occurs very 
near 2 A due to the logarithmic singularity in U'(u>). In 
practice, however, the resonance at weak coupling can 
hardly be observed because the residue of the peak in 
the spin susceptibility Z s is exponentially small. 




CO/CO 2 



FIG. 4. Imaginary part of the dynamical spin susceptibility 
in the superconducting state for T <C A and A = 0.5,1,2, 
determined from the full solution of the Eliashberg equations. 
Dashed line - the normal state result shown for comparison. 
Observe that the resonance peak gets sharper when it moves 
away from 2A. 

Fig.^ shows our results for x(Q,o>) obtained from the 
full solution of the set of three coupled equations at T w 
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and A = 0.5, 1, 2. We clearly see that for A > 1, the spin 
susceptibility has a sharp peak at u — A s . The peak gets 
sharper when it moves away from 2A. At the same time, 
for A = 0.5, which models weak coupling, the peak is very 
weak and is washed out by a small thermal damping. In 
this case, x" omv displays a discontinuity at 2A. 

Before we proceed with the analysis of fermionic prop- 
erties, we show that the spin resonance does not exist 
for s— wave superconductors. In the latter case, the spin 
polarization operator is given by almost the same ex- 
pression as in (p"3|), but with a different sign for <& n <& n+m 
term. The sign difference comes from the fact that the 
two fermions in the spin polarization bubble differ in mo- 
mentum by Q, and the d-wave gap changes sign under 
k + Q. 

One can immediately check using that for a dif- 
ferent sign of the anomalous term, II" does not possess 
a jump at 2A (the singular contributions from £ n I]„ +m 
and cancel each other). Then II' (w) does not 

diverge at 2 A, and hence there is no resonance at weak 
coupling. Still, however, one could naively expect the res- 
onance at strong coupling as at small frequencies II' (lj) 
is quadratic in u> simply by virtue of the existence of the 
threshold for II". It turns out, however, that the reso- 
nance in the case of isotropic s-wave pairing is precluded 
by the fact that II(u; = 0) < (recall that in case of 
d-wave pairing, II(a; = 0) = 0). This negative term over- 
shadows the positive oJ 2 term in II(cj) such that for all 
frequencies below 2A, U(uj) < @. That II (w = 0) < 
in s— wave superconductors can be easily explained. In- 
deed, a negative 11(0) implies that the spin correlation 
length decreases as the system becomes superconducting. 
This is exactly what one should expect as s-wave pair- 
ing involves fermions both from opposite and from the 
same sublattice. The pairing of fermions from the same 
sublattice into a spin-singlet state obviously reduces the 
antiferromagnetic correlation length. 
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FIG. 5. Inelastic neutron scattering intensity for momen- 
tum Q = (n, 7r) as function of frequency for YBa2Cu30e.9. 
Data from Ref. ||. 

The full analysis of the resonance peak requires more 
care as (i) the peak is only observed in two-layer ma- 
terials, and only in the odd channel, (ii) the momen- 
tum dispersion of the peak is more complex than that 
for magnons |3(j], (iii) the peak broadens with under- 
doping jl9 20 1, (iv) in underdoped materials, the peak 



emerges at the onset of the pseudogap, and only sharpens 
at T c |2^,^|. All these features can be explained within 
the spin-fermion model as well. The discussion of these 
effects, however, is beyond the scope of the present pa- 
per. For our present purposes, it is essential that the res- 
onance peak has been observed experimentally, and that 
its frequency for optimally doped materials is around 40 
meV. 

We now proceed to the detailed analysis how the spin 
resonance at A s affects fermionic observables. In each 
instance, we first briefly discuss the results for a d— wave 
gas, and then focus on the strong coupling limit. 



IV. "FINGERPRINTS" OF THE SPIN 
RESONANCE IN FERMIONIC VARIABLES 



C. Comparison with experiments 

In Fig.|^ we show the representative experimental re- 
sult for x"(Q,o;) for optimally doped YB12CU3O6.9 [ j35[ . 
We clearly see a resonance peak at w ~ 41 meV. A 
very similar result has been recently obtained for Bi2212 
pjj . In the last case, the resonance frequency is 43 meV. 
With underdoping, the measured resonance frequency 
goes down |l^j2^]. In strongly underdoped YBi 2 Cu306.6, 
it is approximately 25 meV |L9]]. The very existence of 
the peak and the downturn renormalization of its posi- 
tion with underdoping agree with the theory. Further- 
more, the measured x"(Q> w ) displays an extra feature 
at 60 — 80 meV p0lp5| , which can be possibly explained 
as a 2A effect. 



A. The spectral function 

We first consider the spectral function A^uo) = 
(l/n)\G£(u))\. In the superconducting state, for quasi- 
particles near the Fermi surface 



Ak(uj > 0) = — Im 

7T 



E(6»,w) + e k 



(19) 



where, we remind, T,(6,uj) is the fermionic self-energy 
(which includes a bare oj term), and $(#, u>) is a d— wave 
anomalous vertex. Both can be taken at the Fermi sur- 
face, and generally depend on the direction of kp, which 
we labeled as 9. Also, by definition, Ak(—u>) = A^(oj). 
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1. The d-wave gas 

In a Fermi gas with d-wave pairing, lo) — lo, and 
&(9,w) — A(9) cx cos (29). The spectral function then 
has a 5— functional peak at lo = (A 2 (9) + sl) 1 ^ 2 . It is 
obvious but essential for comparison with the strong cou- 
pling case that the peak disperses with k and far away 
from the Fermi surface recovers the normal state disper- 
sion. 



2. Strong coupling behavior 

Here we consider the spectral function, Ak(u>), for 
fcrmions located near hot spots. As we discussed, for 
those fermions, one can obtain a closed set of equa- 
tions (Eqs. |ll[Jl3| ) which determine A^(uj). Since spin 
fluctuations in a superconducting phase are propagat- 
ing quasiparticles, the effect of the spin scattering on 
fermions near hot spots should be exactly the same as 
the effect of phonon scattering in an s-wave supercon- 
ductor, i.e., in addition to a peak at lo = A, the spectral 
function A^ hs (lo) = A(lo) should possess a singularity at 
lo = ujq = A + A s . 

The behavior of A(lo) near the singularity is very ro- 
bust and can be obtained even without a precise knowl- 
edge of the frequency dependence of S(w) and <&(u>). Our 
reasoning here parallels the one for II(w): all we need to 
know is that near lo = A, T> 2 (lo) — $ 2 (lo) oc lo — A. 
Substituting this form into ( |l3| ) and using the spectral 
representation, we obtain for lo = loq + e 



dx 



o (x(e-z)) 1/2 



(20) 



Evaluating the integral, we find that S" undergoes a fi- 
nite jump at lo = ujq. By Kramers-Kronig relation, this 
jump gives rise to a logarithmic divergence of Exactly 
the same singular behavior holds for the anomalous ver- 
tex &(lu), with exactly the same prefactor in front of the 
logarithm. The last result implies that E(o>) — $(w) is 
non-singular at lo = luq. Substituting these results into 
(|l9|), we find that the spectral function A(lo) emerges at 
lo > loq as l/log 2 (tj — loq), i.e., almost discontinuously. 
Obviously, at a small but finite T, the spectral function 
should have a dip very near lo = lo$, and a hump at a 
somewhat higher frequency. 
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FIG. 6. (a) The quasiparticle spectral function deter- 
mined by solving the coupled Eliashberg equations. The 
peak-dip-hump structure of A(lo) is clearly visible, (b) The 
real (dashed line) and imaginary (solid line) parts of the elec- 
tronic self-energy, (c) The derivative of the imaginary part of 
the self-energy over lo. The singularities at both A + A s and 
at 3A are clearly seen. 

In Fig.|]we present S(w) and A(lo) from a solution of 
the set of three coupled Eliashberg equations at T « 0. 
This solution is consistent with our analytical estimate. 
We clearly see that the fermionic spectral function has 
a peak-dip-hump structure, and the peak-dip distance 
exactly equals A s . We also see in Fig.||that the spectral 
function is nonanalytic at lo = 3A. As we discussed, 
this nonanalyticity is peculiar to the spin-fermion model, 
and is due to the nonanalyticity of the dynamical spin 
susceptibility at lo = 2 A. 

Another "fingerprint" of the spin-fluctuation scatter- 
ing can be observed by studying the evolution of the spec- 
tral function as one moves away from the Fermi surface. 
The argument here goes as follows: at strong coupling, 
where A > lv s {, probing the fermionic spectral function 
at frequencies progressively larger than A, one eventually 
probes the normal state fermionic self-energy atw> Lo s f. 
Due to strong spin-fermion interaction, this self-energy 
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is large. Indeed, the solution of Eq.(13) with <f> = and 
H(w) = iu yields at T = ||§ 



maximum at o> = £^/oj. The positions of the peak and 
the dip versus k — kp are presented in FigMj 



EH=«(1+ ^-j ) (21) 

V i + (i-^) 1/2 / 

For 3> w s /j E(tt>) becomes 

S(tj) = cj(l + (^/|t i ;|) 1 / 2 ) (22) 

Substituting this form into the fermionic propagator, we 
find that up to u> ~ Q, the spectral function in the nor- 
mal state does not have a quasiparticle peak at ui = £k- 
Instead, it only displays a broad maximum at u> = 
Alternatively speaking, at uj s f < lu < Co, the spectral 
function in the normal state displays a non- Fermi liquid 
behavior with no quasiparticle peak (see Fig. This 
particular non-Fermi liquid behavior (S(w) oc y/uj) is as- 
sociated with the closeness to a quantum phase transition 
into an antiferromagnetic state. Indeed, at the transition 
point, uj s f = 0, and hence £(oj) cx -y/uJ extends down to a 
zero frequency. 

The absence of the quasiparticle peak in the normal 
state implies that the sharp quasiparticle peak which we 
found at u> — A for momenta at the Fermi, cannot simply 
disperse with k, as it does for noninteracting fcrmions 
with a ci-wave gap. Specifically, the quasiparticle peak 
cannot move further in energy than A + A s as at larger 
frequencies, the spin scattering becomes possible, and the 
fermionic spectral function should display roughly the 
same non-Fermi-liquid behavior as in the normal state. 
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FIG. 7. The normal state spectral function, relevant for the 
high frequency behavior in the superconducting state. Note 
the absence of a quasiparticle peak. This is the consequence of 
the proximity to an antiferromagnetic quantum critical point. 

In Fig.||a we present successive plots for the spectral 
function as the momentum moves away from the Fermi 
surface. We see exactly the behavior we just described: 
the quasiparticle peak cannot move further than A + 
A s . Instead, when k — kp increases, it gets pinned at 
A + A s and gradually looses its spectral weight. At the 
same time, the hump disperses with k and for frequencies 
larger than A + A s gradually transforms into a broad 
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FIG. 8. a) Frequency dependence of the spectral function 
in the superconducting state for different e^. The curve at 
the bottom has a highest e^. No coherent quasiparticle peak 
occurs for energies larger than A + A 3 . Instead, the spectral 
function displays a broad maximum, similar to that in the 
normal state. 

3. Comparison with the experiments 

The quasiparticle spectral function at various mo- 
menta is measured in angle resolved photoemission 
(ARPES) experiments. In a sudden approximation (an 
electron, hit by light, leaves the crystal without fur- 
ther interactions with other electrons and without pay- 
ing attention to selection rules for the optical transition 
to its final state), the photoemission intensity is given 
by I\t{ui) = Ak(Lu)np(uj) where rip is the Fermi func- 
tion. ARPES data of Ref. § for near optimally doped, 
T c = 87K in Bi2212 for momenta near a hot spot are pre- 
sented in Fig|]. We clearly see that the intensity displays 
a peak/dip/hump structure. A sharp peak is located at 
~ 40meV, and the dip is at 80meV. 
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FIG. 9. ARPES spectrum for near optimally doped Bi2212 
for momenta close to the hot spots. Data from Ref. [6]. 
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The experimental peak-dip distance is 42meV . The 
neutron scattering data plj] on Bi2212 with nearly the 
same T c = 91K yielded A s = 43meV which is in excellent 
agreement with the ARPES data. Furthermore, with un- 
derdoping, the peak-dip energy difference decreases and, 
up to error bars, remains equal to A s . This behavior is 
illustrated in Fig.[lC|. 
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FIG. 10. The experimental peak-dip distance at various 
doping concentrations vs A a extracted from neutron measure- 
ments. Data from Ref.[30] 

Finally, in Fig.[H] we present experimental results for 
the variation of the peak and hump positions with the 
deviation from the Fermi surface. We clearly see that the 
hump disperses with k — kp and eventually recovers the 
position of the broad maximum in the normal state. At 
the same time, the peak shows little dispersion, and does 
not move further in energy than A + A s . Instead, the 
amplitude of the peak just dies off as k moves away from 
kp. This behavior is fully consistent with the theory. 
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FIG. 11. The experimental peak and hump positions with 
the deviation from the Fermi surface. We clearly see that 
the hump disperses with k — kF and eventually recovers the 
position of the broad maximum in the normal state, while the 
peak position changes little with the deviation from kF. Data 
from Ref.[6]. 

We regard the presence of the dip at A + A s , and the 
absence of the dispersion of the quasiparticle peak are two 
major "fingerprints" of strong spin-fluctuation scattering 
in the spectral density of cuprate superconductors. 



B. The density of states 

The quasiparticle density of states, N(ui), is the mo- 
mentum integral of the spectral function: 



N(uj) 



4tt 



2 A k (w). 



(23) 



Substituting Ak(uj) from Eq.(|19|) and integrating over £k, 
one obtains 



2tt 



where, we remind, is the angle along the Fermi surface. 
As before, we first consider N(oj) in a ei-wave gas, and 
then discuss strong coupling effects. 



1. Density of states in a d-wave gas 

Consider for simplicity a circular Fermi surface for 
which Ak = A cos (26). Integrating in (|24|) over we 
obtain lEl 



N(uj) = Rc 



uj 
2^ 



dd 



y/u 2 - A 2 cos 2 (26>) 

2 f K(A/u) for uj > A 

7T 1 (lu/A)K{uj/A) for u < A ' 



(25) 



where K(x) is the elliptic integral of first kind. We see 
that N(ui) <~ ui for uj <S A and diverges logarithmically as 
(1/tt) ln(8A/| A — ui\) for uj ps A. At larger frequencies, 
N(uj) gradually decreases to a frequency independent, 
normal state value of the DOS, which we normalized to 
1. The plot of N(uj) in a d— wave BCS superconductor is 
presented in Fig 112] 
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FIG. 12. Density of states of a noninteracting Fermi gas 
with d-wave gap (solid line) and with s-wave gap (dashed 
line) . 

For comparison, in an s-wave superconductor, the DOS 
vanishes at lo < A and diverges as {uj — A) -1 / 2 at lo > A. 
We see that a d-wave superconductor is different in that 
(i) the DOS is finite down to the smallest frequencies, 
and (ii) the singularity at uj = A is weaker (logarithmic) . 
Still, however, N{lo) is singular only at a frequency which 
equals to the largest value of the d— wave gap. This il- 
lustrates a point we made earlier that the angular de- 
pendence of the d-wave gap reduces the strength of the 
singularity, but does not wash it out over a finite fre- 
quency range. 



2. Density of states at strong coupling 

We first show that the linear behavior of N{lo) at low 
frequencies and the logarithmic divergence at lo = A, ob- 
served in a gas, are in fact quite general and are present 
in an arbitrary d— wave superconductor. Indeed, at low 
frequencies ReS {8, lo) oc uj, Re<fr{8,co) oc (9— 6> no dc) where 
#nodc = f ) T"- are the positions of the node of the d- 



wave gap. Substituting these forms into (|24|) and inte- 
grating over 9 we obtain N(uj) oc lo. Similarly, expanding 
5]2 _ <j)2 near a S p t ; where, at least at strong cou- 
pling, the d— wave gap is at maximum | p^ j, we obtain 
Y?{8, uj) - $ 2 (0, lo) cx {uj - A) + B8 2 , where 9 = 8- 9 hs , 
and B > 0. Here, 6>h s is the position of a hot spot on 
the Fermi surface. The integration over 9 then yields 



N{lo) oc Re 



d8 



ioRin-Ai 



B8 2 + (fi — A) 



(26) 



This result implies that at strong coupling, the DOS in a 
d— wave superconductor still has a logarithmic singularity 
at uj = A, although the overall factor for the singular 
piece depends on the strength of the interaction. 

We see therefore that the behavior of the density of 
states in a d-wave superconductors for lo < A is quite 
robust against strong coupling phenomena. This nice 
universality on the other hand makes N{lo < A) not 
very sensitive to a specific mechanism for the pairing. 



We now demonstrate that at strong coupling, the DOS 
possesses extra peak-dip features, associated with the sin- 
gularities in S(w) and $(w) at lo = uj = A + A s . The 
analytical consideration proceeds as follows Con- 
sider first a case when the gap is totally flat near a hot 
spot. At co = ujo, both £(w) and $(w) diverge logarithmi- 
cally. Substituting these forms into (|24|), we immediately 
obtain that N{lo) has a logarithmic singularity: 



N sing {uj) oc log 



1 



\UJ - LO \ 



1/2 



(27) 



This singularity gives rise to a strong divergence of 
dN{uo)/du) at uj = loq- This behavior is schematically 
shown in Fig. (^)a. In part (b) of this figure we 
present the result for N{uo) obtained by the solution of 
Eliashberg-type Eqs. |ll"||T3| . A small but finite tempera- 
ture was used to smear out divergences. We recall that 
the Eliashberg set does not include the angular depen- 
dence of the gap near hot spots, and hence the result for 
the DOS should be compared with Fig. (|l3|)a. We clearly 
see that N{uj) has a second peak at uj = loq. This peak 
strongly affects the frequency derivative of N{lo) which 
has a predicted singular behavior near loq. 

Note also that a relatively small magnitude of the sin- 
gularity in N{lo) is a consequence of the linearization of 
the fermionic dispersion near the Fermi surface. For ac- 
tual efc chosen to fit ARPES data , the nonlinearities in 
the fermionic dispersion occur at energies comparable to 
loo. This is due to the fact that hot spots are located close 
to (0, 7r) and related points at which the Fermi velocity 
vanishes. As a consequence, the momentum integration 
of the spectral function should have less drastic smearing 
effect than in our calculations, and the frequency depen- 
dence of N{lo) should more resemble that of A{lo) for 
momenta where the gap is at maximum. 
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FIG. 13. (a) The behavior of the SIN tunneling conduc- 
tance (i.e., DOS) in strongly coupled d-wave superconductor. 
Main pictures - N(u>), insets - dN{uj)/du. (a) The schematic 
behavior of the DOS for a flat gap. (b) The solution of the the 
set of the Eliashberg-type equations for flat gap. The shaded 
regions are the ones in which the flat gap approximation is 
incorrect as the physics is dominated by nodal quasiparticles. 
(c) The schematic behavior of N(u>) for the quadratic varia- 
tion of the gap near its maxima, (d) The expected behavior 
of the DOS in a real situation when singularities are softened 
out by finite T or impurity scattering. The position of A + A s 
roughly corresponds to the minimum of dN(tv)/dw. 



For real materials, in which singularities are removed 
by e.g., impurity scattering, N(uj) likely has a dip at 
uj > luq and a hump at a larger frequency. This is 
schematically shown in Fig. (13)d. The location of ujq is 
best described as a point where the frequency derivative 
of the DOS passes through a minimum. 

The singularity in dN(uj)/du at wq, and the dip- hump 
structure of N(uj) at u> > uj are another "fingerprints" 
of spin-fluctuation mechanism in the single particle re- 
sponse. 



For momentum dependent gap, the behavior of 
fermions near hot spots is the same as when the gap 
is flat, but now ojq depends on 9 as both A and A s vary 
as one moves away from a hot spot. The variation of 
A is obvious, the variation of A s is due to the fact that 
this frequency scales as A 1 / 2 . Since both A and A s are 
maximal at or very near a hot spot, we can write quite 
generally 



w 



wo — a9 



(28) 



where a > 0. The singular pieces of the self-energy and 
the anomalous vertex then behave as log(o;o— uj — a6 ) . 
Substituting these forms into ( |2~4| ) and using the fact that 
— $(<a>) ~ const at w k luq, we obtain 



N s i ng {uj) oc Re 



dO log(w 



uj - 



a6 2 y 



-1/2 



(29) 



A straightforward analysis of the integral shows that now 
N(uj) has a one-sided nonanalyticity at uj = luq: 



-WsingM = -BQ(u - UJ ) 



UJ — UJq 



log(w - Uj a )\ 



1/2 



(30) 



where B > 0, and 0(x) = 1 for x > 0, and Q(x) = for 
x < 0. This nonanalyticity gives rise to a cusp in N(uj) 
right above ujq, and a one-sided square-root divergence 
of the frequency derivative of the DOS. This behavior 
is shown schematically in Fig. fll3|) c. Comparing this 
behavior with the one in Fig. (|l3[)a for a flat gap, we 
observe that the angular dependence of the gap predom- 
inantly affects the form of N(u>) at uj < ujq. At these 
frequencies, the angular variation of the gap completely 
eliminates the singularity in N(lu). At the same time, 
above ujq, the angular dependence of the gap softens the 
singularity, but still, the DOS sharply drops above ujq 
such that the derivative of the DOS diverges at approach- 
ing ujq from above. Alternatively speaking, in a d— wave 
superconductor, the singularity in the DOS is softened by 
the angular dependence of the gap, but it still holds at 
a particular frequency related to the maximum value of 
the gap. This point is essential as it enables us to read off 
the maximum gap value directly from the experimental 
data without "deconvolution" of momentum averages. 



3. Comparison with the experiments 

As we already mentioned, the fermionic DOS N(uj) 
is proportional to the dynamical conductance dl/dV 
through a supcrconductor-insulator-normal metal (SIN) 
at uj = eV [Q. The drop in the DOS at ujq can be refor- 
mulated in terms of SIN conductance as follows: if the 
voltage for SIN tunneling is such that eV — ujq, then an 
electron which tunnels from the normal metal, can emit 
a spin excitation and fall to the bottom of the band loos- 
ing its group velocity. This obviously leads to a sharp 
reduction of the current and produces a drop in dl/dV. 
This process is shown schematically in Figjl4| ]To| . 
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FIG. 14. The schematic diagram for the dip features in 
the SIN and SIS tunneling conductances (figures a and b, re- 
spectively) . For the SIN tunneling, which measures fermionic 
DOS, the electron which tunnels from a normal metal can 
emit a propagating magnons if the voltage eV = A + A s . 
After emitting a magnon, the electron falls to the bottom of 
the band. This leads to a sharp reduction of the current and 
produces a drop in dl/dV. For SIS tunneling, the physics 
is similar, but one first has to break an electron pair, which 
costs energy 2A. 

The SIN tunneling experiments have been performed 
on YBCO and Bi2212 materials We reproduce these 
data in Fig. |l5|. Very similar results have been recently 
obtained by Pan et al [ pT| . 
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FIG. 15. The experimental result for the differential con- 
ductance through the SIN tunneling junction. This conduc- 
tance is proportional to the quasiparticle DOS. The data are 
for Bi2212 and are taken from Ref. Q. 

At low and moderate frequencies, the SIN conduc- 
tance displays a behavior which is generally expected in 
a d— wave superconductor, i.e., it is linear in voltage for 
small voltages, and has a peak at eV = A where A is 
the maximum value of the d— wave gap [jl4|,[ll). The 
value of A extracted from tunneling agrees well with 
the maximum value of the gap extracted from ARPES 
measurements JT^Jl3|]. At frequencies larger than A, the 
measured SIN conductance clearly displays an extra dip- 
hump feature which become visible at around optimal 
doping, and grows in amplitude with underdoping [flif . 
At optimal doping, the distance between the peak at A 
and the dip is around 40meV. This is consistent with 
A s extracted from neutron measurements. The doping 
dependence of the peak-dip distance, and the behavior of 
dN(ui)/dui have not been studied in detail, to the best of 
our knowledge. This analysis is clearly called for. 



C. SIS tunneling 

The measurements of the dynamical conductance 
dl/dV through a superconductor - insulator - supercon- 
ductor (SIS) junction is another tool to search for the 
fingerprints of the spin-fluctuation mechanism. The con- 
ductance through this junction is the derivative over volt- 
age of the convolution of the two DOS Q: dl/dV oc S(lo) 
where 



dnN(iu - n) d n N(n) 



(31) 



1. SIS tunneling in a d-wave gas 

The DOS in a d— wave gas is given in Eq. (p5|) . Sub- 
stituting this form into Eq. ( |3l| ) and integrating over fre- 
quency we obtain the result presented in Fig.|l6|. At small 
lu, S(u>) is quadratic in frequency p4j. This is an obvious 
consequence of the fact that the DOS is linear in w. At 
lu = 2A, S(lo) undergoes a finite jump. This jump is re- 
lated to the fact that near 2 A, the integral over the two 
DOS includes the region around SI = A where both JV(f2) 
and N((jj — Q) are logarithmically singular, and dnN(£i) 
diverges as l/(f2 — A) . The singular contribution to S(w) 
from this region can be evaluated analytically and yields 



1 



S(u>) = --»P 



dx log | £ | 1 

7nr = — TrSignfw — 2A) 

x + uj-2A 2 6 v ; 



(32) 



Observe that the amplitude of the jump in the SIS con- 
ductance is a universal number which does not depend 
on the value of A. 




FIG. 16. The SIS tunneling conductance, dl/dV, in a 
d— wave BCS superconductor. The dashed line shows 
for s-wave case for comparison. 

At larger frequencies, S(u>) continuously goes down 
and eventually approaches a value of S(uj — * oo) = 1. 



2. SIS tunneling at strong coupling 

As in the previous subsection, we first demonstrate 
that the quadratic behavior at low frequencies and the 
discontinuity at 2A survive at arbitrary coupling. In- 
deed, the quadratic behavior at low u> is just a conse- 
quence of the linearity of N(u) at low frequencies. As 
shown above, this linearity is a general property of a 
d— wave superconductor. Similarly, the logarithmic di- 
vergence of the DOS at u> — A causes the discontinuity in 
the SIS conductance by the same reasons as in a d— wave 
gas. 

We next consider how the singularity in S(w) at loq 
affects S(u>). From a physical perspective, we should ex- 
pect a singularity in S(u>) at u> — A + ujq = 2 A + A s . 



13 



Indeed, to get a nonzero SIS conductance, one has to 
first break a Cooper pair, which costs an energy of 2A. 
After a pair is broken, one of the electrons becomes a 
quasiparticle in a superconductor and takes the energy 
A, while the other tunnels. If eV = A + ujq, the electron 
which tunnels through a barrier has energy u>o, and can 
emit a spin excitation and fall to the bottom of the band 
(see Fig. |l4|). This should produce a drop in dl/dV by 
the same reasons as for SIN tunneling. This behavior is 
schematically shown in Fig. p"7| . 

Consider this effect in more detail. We first notice 
that lj = A + ujq is special for Eq.(|3l|) because both 
dN(Cl)/dfl and N(oj — fl) diverge at the same energy, 
O = luq. Substituting the general forms of N(u>) near 
lu = ljq and to = A, we obtain after simple manipulations 
that for a flat gap, S(uj) has a one-sided divergence at 
oj = uj + A = 2A + A s . 



5 (e)oc 



9( 



(33) 



where e = lu — (luq + A). This obviously causes the 
divergence of the frequency derivative of S(uj) (i.e., of 
d 2 I/dV 2 ). This behavior is schematically shown in Fig. 
p7|a. In Fig. [T^ b we present the results for S(uj) ob- 
tained by integrating theoretical N(lu) from the previous 
subsection (see Fig.|I3|b). We clearly see that S(lu) and 
its frequency derivative are singular at lu = 2 A + A s , in 
agreement with the analytical prediction. 

For quadratic variation of the gap near the maxima, 
the calculations similar to those for the SIN tunneling 
yield that S(u>) is continuous through 2A + A S) but its 
frequency derivative diverges as 



dSjfj) 

dbJ 



(X P 



dx 



o (x\logx\y/ 2 (x-e) \elog\e\\^- 



(34) 



The singularity in the derivative implies that near e = 



S{e) = 5(0) - C 9(-e) 



log(- 



1/2 



(35) 



where C > 0. This behavior is schematically presented 
in Fig. fl7|d. We again see that the angular dependence 
of the gap softens the strength of the singularity, but the 
singularity remains confined to a single frequency lu — 
2A + A s . 

In real materials, the singularity in S(lu) is softened 
and transforms into a dip slightly below 2 A + A s , and 
a hump at a frequency larger than 2 A + A s . The fre- 
quency 2A + A s roughly corresponds to a maximum of 
the frequency derivative of the SIS conductance. 
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FIG. 17. (a) The schematic behavior of the SIS tunneling 
conductance, S(u), in a strongly coupled d-wave supercon- 
ductor. Main pictures - S(uj), insets - dS(ui) / 'duo. (a) The 
schematic behavior of S{uj) for a flat gap. (b) The solution of 
the the set of the Eliashberg-type equations for flat gap using 
the DOS from the previous subsection. The shaded regions 
are the ones in which the flat gap approximation is incorrect 
as the physics is dominated by nodal quasiparticles. (c) The 
schematic behavior of S(oj) for quadratic variation of the gap 
near its maxima, (d) The expected behavior of the SIS con- 
ductance in a real situation when singularities are softened 
out by finite T or by impurity scattering. 2A + A 3 roughly 
corresponds to the maximum of dS(oj) / du. 



3. comparison with experiments 

Recently, Zasadzinski et al. |l5| obtained and carefully 
examined their SIS tunneling data for a set of Bi2212 
materials ranging from overdoped to underdoped |L5| |. 
Their data, presented in Fig.Q clearly show that be- 
sides a peak at 2 A, the SIS conductance also has a dip 
and a hump at larger frequencies. The distance between 
the peak and the dip (« A s in our theory) is close to 2A 
in overdoped Bi2212 materials, but goes down with un- 
derdoping. Near optimal doping, this distance is around 
40mcV. For underdoped, T c = 74K material, the peak- 
dip distance is reduced to about 30meV. 
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FIG. 18. SIS tunneling data from Ref.[ll] for a set of 
Bi2212 materials ranging from overdoped to underdoped. 

These results are in qualitative and quantitative agree- 
ment with ARPES and neutron scattering data, as well 
as with our theoretical estimates. The most important 
aspect is that with underdoping, the experimentally mea- 
sured peak-dip distance progressively shifts down from 
2 A. This is the key feature of the spin-fluctuation mech- 
anism, We regard the experimental verification of this 
feature in the SIS tunneling data is a strong argument in 
favor of a magnetic scenario for superconductivity. 



D. Optical conductivity and Raman response 

Another observables sensitive to ujq are the optical con- 
ductivity and the Raman response. Both are propor- 
tional to the fully renormalized particle-hole polarization 
bubble, but with different signs attributed to the bubble 
made of anomalous propagators. Namely, after integrat- 
ing in the particle- hole bubble over £k, one obtains 



R(uj) =Im / duj'd9V 2 (9)U r (9,uj,uj') 

o-(w)=Im — l — [ duj'd9n a (9,uj,u') 
uj + id J 



(36) 
(37) 



where V(0) is a Raman vertex which depends on the 
scattering geometry Q, and 



U r ,a(Q, u, uj') 



£>+£>_(£>+ + £>_) 



(38) 



Here a = — 1 for II r , and a = 1 for H a . Also, E± 
V (uj±) and $± = $ (oj±) with uj± = uj' ± uj/2 as well as 
D± = ($| - El) 1 / 2 . Note, £ and $ depend on uj and 9. 



1. Optical and Raman response in a d-wave gas 

In a superconducting gas, the optical conductivity van- 
ishes identically for any nonzero frequency due to the 
absence of a physical scattering between quasiparticles 
in a gas. The presence of a superconducting condensate, 
however, gives rise to a 5 functional term in a at uj = 0: 
<j(uj) = ttS(uj) J d9duj'H a (9,0 7 uj'). This behavior is typi- 



cal for any BCS superconductor 39 1, and holds for both 
s— wave and d— wave superconductors. The behavior of 
<j(uj) in a superconducting gas with impurities, causing 
inelastic scattering, is more complex and has been dis- 
cussed for a d— wave case in, e.g., Ref. plL 

The form of the Raman intensity depends on the scat- 
tering geometry. For mostly studied B\ g scattering, the 
Raman vertex has the same angular dependence as the 
d-wave gap, i.e., V(9) 



40 1 . Straightfor- 



oc cos (26) 

ward computations then show that at low frequencies, 
R(lu) cx uj 3 [^0). For a constant V(9), we would have 
R(uj) cx uj. 

Near u> = 2A, the B\ g Raman intensity is singular. For 
this frequency, both D + and D_ vanish at uj' = and 
9 = 0. This causes the integral for R(uj) to be divergent. 
The singular contribution to R(uj) can be obtained ana- 
lytically by expanding in the integrand to leading order 
in uj' and in 6. Using the spectral representation, we then 
obtain, for uj = 2A + 5 ft 



R(uj) = dn d6 



1 



1 



(39) 



Where, as before, 8 = 9 — 9h s For a flat band (a = 0), 
R(uj) oc Re[(uj — 2A) -1 / 2 ]. For a^O, i.e., for a quadratic 
variation of the gap near its maximum, the 2D integra- 
tion in (|39|) yields R(uj) cx | log e | . At larger frequencies 
R{uj) gradually decreases. 

The behavior of R(uj) in a gas is shown in Figjl9|. 
Observe that due to interplay of numerical factors, the 
logarithmic singularity shows up only in a near vicinity 
of 2 A, while at somewhat larger uj, the angular depen- 
dence of the gap becomes irrelevant, and R(uj) behaves 
as (uj — 2A) 1 / 2 , i.e., as for a flat gap. 
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FIG. 19. The behavior of the Raman response in a BCS su- 
perconductor with a flat gap (dashed line), and for a quadratic 
variation of the gap near its maximum (solid line). 

2. Raman and optical response at strong coupling 

A nonzero fermionic self-energy mostly affects the op- 
tical conductivity for a simple reason that it becomes 
finite in the presence of the spin scattering which can re- 
lax fermionic momentum. For a momentum-independent 
gap, a finite conductivity emerges above a sharp thresh- 
old. This threshold stems from the fact that at least one 
of the two fermions in the conductivity bubble should 
have a finite £", i.e., its energy should be larger than 
loq. Another fermion should be able to propagate, i.e., 
its energy should be larger than A. The combina- 
tion of the two requirements yields the threshold for 
a(u> > 0) at 2A + A s , i.e., at the same frequency where 
the SIS tunneling conductance is singular. One can easily 
demonstrate that for a flat gap, the conductivity emerges 
above the threshold as e 1 / 2 /log 2 e, where, we remind, 
e = lu — A — luq = uj — (2 A + A s ). This singularity ob- 
viously causes a divergence of the first derivative of the 
conductivity at e = +0. 

For a true d— wave gap, the conductivity is finite for 
all frequencies simply because the angular integration in 
Eq.([37j) involves the region near the nodes, where £" is 
nonzero down to the lowest frequencies. Still, however, 
we argue that the conductivity is singular at luq + A. 
Indeed, replacing, as before, e by e + 2 , substituting this 
into the result for the conductivity for a flat gap, and 
integrating over 9, we find that for momentum dependent 
gap, the conductivity itself and its first derivative are 
continuous at e = 0, but the second derivative of the 
conductivity diverges as d 2 a/duj 2 cx l/(|e| log 2 e). 

In Fig|2(] we show the result for conductivity obtained 
by solving the set of coupled Eliashberg-type equations 
Eqs. |ll[jl3| . We clearly see an expected singularity at 
2 A + A s . Note by passing that at higher frequencies, 
the theoretical is inversely linear in u> 0). 
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FIG. 20. Frequency dependence of the optical conductivity 
computed from the self energy and pairing vertex determined 
from the Eliashberg equations valid near hot spots. The onset 
of the optical response is lo — 2A + A s . The contribution from 
the nodes (not included in calculations) yields a nonzero con- 
ductivity at all u) and softens the singularity at uj — 2 A + A s . 
The insert shows the experimental data of Puchkov et al. jl6| . 

For the Raman intensity, the strong coupling effects are 
less relevant. First, one can prove along the same lines as 
in previous subsections that the cubic behavior at low fre- 
quencies for Big scattering (and the linear behavior for 
angular independent vertices), and the logarithmic sin- 
gularity at 2A are general properties of a d— wave super- 
conductor, which survive for all couplings. Thus, similar 
to the density of states and the SIS-tunneling spectrum, 
the Raman response below 2A is not sensitive to strong 
coupling effects. Second, near ujq + A, singular contribu- 
tions which come from £+£_ and $ + $_ terms in II r in 
Eq.(^) cancel each other. As a result, we found that for 
a flat gap, only the second derivative of R(u>) diverges 
at A + luq. For a quadratic variation of a gap near its 
maximum, the singularity is even weaker and shows up 
only in the third derivative of R(oj). Obviously, this is 
a very weak effect, and its determination requires a high 
quality of the experiment. Notice, however, that, as we 
already mentioned, due to the closeness of hot spots to 
(0, 7r) and related points, where vp vanishes, the smearing 
of the singularity due to momentum integration may be 
less drastic than in our theory where we used a linearized 
fermionic dispersion with some finite Fermi velocity. 

3. Comparison with experiments 

Our theoretical considerations show that optical mea- 
surements are much better suited to search for the "fin- 
gerprints" of a magnetic scenario, then Raman measure- 
ments. Evidences for strong coupling effects in the opti- 
cal conductivity in superconducting cuprates have been 
reported in Refs. |l^,[ll||^|,[l6| . We present the experi- 
mental data for <j{u>) in optimally doped YBCO in the 
inset of Fig.|2(]. We see that the conductivity drops at 
about lOOmeV. Earlier tunneling measurements of the 
gap in optimally doped YBCO yielded A = 29meV [0. 
Combining this with A s = 41meV, we find 2A + A s 
lOOmeV, consistent with the data. We consider this 
agreement as another argument in favor of a magnetic 
scenario. 
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FIG. 21. A calculated frequency dependence ol 

W{uj) — ^j(wRe(j~ 1 (cj)). This quantity is a sensitive mea- 
sure for fine structures in the optical response. 
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for = 

of the deep minimum agrees well with 2 A + A s . The extrema 
at higher frequencies are consistent with 4A and 2 (A + A 
predicted by the theory. 



4- Fine structure of optical conductivity 

We now argue that the measurements of optical con- 
ductivity allow one not only to verify the magnetic sce- 
nario, but also to independently determine both A s and 
A in the same experiment. We discussed several times 
above that in a magnetic scenario, the fermionic self- 
energy is singular at two frequencies: at ujq = A + A s , 
which is the onset frequency for spin-fluctuation scat- 
tering near hot spots, and at u> = 3 A, where fermionic 
damping near hot spots first emerges due to a direct 
four-fermion interaction. We argued that in the spin- 
fluctuation mechanism, both singularities are due to the 
same underlying interaction, and their relative intensity 
can be obtained within a model. 

In general, the singularity at 3A is much weaker at 
strong coupling, and can be detected only in the analysis 
of the derivatives of the fermionic self-energy. We remind 
that the singularity in at ojq gives rise to singularity 



in the conductivity at A + luq , while the 3 A singularity 
in obviously causes a singularity in conductivity at 
lu = 4A. Besides, we should also expect a singularity in 
ct(w) at 2ujo, as at this frequency both fermions in the 
bubble have a singular ~E(ujq). 

For phonon superconductors, a fine structure of opti- 
cal conductivity has been analyzed by studying a second 
derivative of conductivity via W(uj) = -^(ujRea" 1 (uj)) 
which is proportional to a 2 (uj)F(uj) where a(u) is an ef- 
fective electron-phonon coupling, and F (u>) is a phonon 

dos a 



In Fig. 21 we present the theoretical result for W(u>) in 
our model |4^|. First, we clearly see that there is a sharp 
maximum in W(u>) near 2A + A s , which is followed by a 
deep minimum. This form is consistent with our analyt- 
ical observation that for a flat gap (which we used in our 
numerical analysis), the first derivative of conductivity 
diverges at u = 2 A + A s . At a finite T which is a nec- 
essary attribute of a numerical solution, the singularity 
is smoothened, and the divergence is transformed into a 
maximum. Accordingly, the second derivative of the con- 
ductivity should have a maximum and a minimum near 
2A + A s . We found from our numerical analysis that the 
maximum shifts to lower frequencies with increasing T, 
but the minimum moves very little from 2 A + A s , and is 
therefore a good measure of a magnetic "fingerprint" . 

Second, we see from Fig.^T] that besides the maximum 
and the minimum near 2A + A s , W{uj) has extra extrema 
at 4 A and 2luq — 2 A + 2A S . These are precisely the 
extra features that we expect, respectively, as a primary 
effect due to a singularity in at oj — 3 A, and as a 

secondary effect due to a singularity in E(w) at w = wq. 

The experimental result for W(u) shown in Fig|2^. We 
see that the theoretical and experimental plots of W(u>) 
look rather similar. Furthermore, the relative intensi- 
ties of the peaks are at least qualitatively consistent with 
the theory. Identifying the extra extrema in the exper- 
imental W(uj) with 4A and 2A + 2A S , respectively, we 
obtain 4A - 130meV, and 2A + 2A S ~ 150meV. This 
yields A ~ 32meV, in good agreement with earlier mea- 
surements [^9|, and A s ~ 45meV, which is only slightly 
larger than the resonance frequency extracted from neu- 
tron measurements [ fl9l]2l| . Despite the fact that the de- 
termination of a second derivative of a measured quantity 
is a very subtle procedure, the very presence of the extra 
peaks and the fact that their positions fully agree with 
the theory, is an indication that the fine structures in 
the optical response may indeed be due to strong spin 
fermion scattering. 



V. COMPARISON WITH OTHER WORKS 

We now discuss how our work is related to other stud- 
ies. As we stated in the very beginning of the paper, 
the fact that the interaction with a bosonic mode with 
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frequency loq gives rise to a fermionic damping in a su- 
perconductor above u>o + A, is known for conventional 
s— wave superconductors. [|[|4) (see also and Q). 
Recent angle-integrated photoemission data for lead and 
niobium |L7| clearly demonstrated that the photoemis- 
sion intensity at low T possesses peak-dip features, but 
the dip is located well above 3 A. 

The reduction of the spin damping below 2A in a 
rf— wave superconductor has been discussed in It 
has also been argued earlier that the interaction with a 
nearly resonant collective mode peaked at Q explains the 
ARPES data. Qualitative arguments for this have been 
displayed by Shen and Schrieffer |5l|| and Norman and 
Ding |fi^| . Norman and Ding also conjectured that the 
peak-dip separation may be related to the frequency of a 
neutron peak, and presented weak-coupling calculations 
for a model in which fcrmions interact with a resonance 
bosonic mode. The results of this analysis agree with the 
ARPES data. From this perspective, the novelty of our 
approach is in that it presents controlled strong-coupling 
calculations for a low-energy spin-fcrmion model, which 
verify and extend earlier ideas. 

It has been also realized earlier that in a d— wave BCS 
superconductor, the dynamical spin response at wave 
vector Q contains an excitonic pole below 2A. This 
has been demonstrated by a number of researchers fl53| . 
The earlier studies, however, considered a weak coupling 
limit, when a bare particle-hole bubble (a building block 
of RPA series) is made out of free fermions. However, as 
we discussed in Sec. [II, performing weak coupling cal- 
culations consistently, one finds A s exponentially close 
to 2A, and the exponentially small residue of the reso- 
nance peak. Our work extends the idea that the neutron 
resonance is a spin exciton to the strong coupling limit. 

The results for the density of states, the SIS tunnel- 
ing conductance, the Raman intensity, and the optical 
conductivity in a d— wave gas have been studied in de- 
tail by a number of researchers H]. In our approach, 
we used these results as a zero-order theory, and consid- 
ered strong coupling feedback effects on top of it. As 
far as we know, for DOS, SIS tunneling and Raman in- 
tensity, these feedback effects have not been studied be- 
fore. The form of the optical conductivity and of W(lS) 
in cuprates has been recently studied by Carbotte and 
co-workers |lj,^,^] . They also argued that the analysis 
of W(lo) supports a magnetic scenario for the pairing. 
There is, however, some discrepancy between Ref. fl8| l 
and out work: it was argued in Jl8| that the largest, 
broad maximum in W(tu) is shifted down from 2A + A s 
due to the angular dependence of the gap, and is located 
at A + A s . We also found that the broad maximum in 
W(uj) is located at a frequency smaller than 2A+ A s . We 
however, attribute this reduction to finite T. We showed 
in the text that for T <C A, the singularity in W(cu) 
in a d— wave superconductor is still located precisely at 
2A + A s . 



Finally, several groups |57j] recently suggested phe- 
nomenologically that there is a connection between neu- 
tron resonance and specific heat anomaly in cuprates [^8) . 
The verification of this connection within the spin- 
fcrmion model is called for. 



VI. CONCLUSIONS 

In this paper we demonstrated that the same feedback 
effect which in the past allowed one to verify the phononic 
mechanism of a conventional, s-wave superconductivity, 
is also applicable to cuprates, and may allow one to 
experimentally detect the "fingerprints" of the pairing 
mechanism in high T c superconductors. We argued that 
for spin-mediated fermionic scattering, and for the hole- 
like normal state Fermi surface, the fermionic spectral 
function, the density of states, the SIS tunneling con- 
ductance, and the optical conductivity are affected in a 
certain way by the interaction with collective spin excita- 
tions which in the superconducting state are propagating, 
magnon-like quasiparticles with the gap A s . We have 
shown that the interaction with propagating spin excita- 
tions gives rise to singularities at frequencies A + A s for 
the spectral function near hot spots and the DOS, and at 
2A + A s for the SIS tunneling conductance and the opti- 
cal conductivity. We demonstrated that the value of A s 
extracted from these experiments agrees well with the re- 
sults of neutron experiments which measure A s directly. 

We further argued that in optical measurements one 
can also detect subleading singularities at 4A and 2A + 
2A S , and that these fine features have been observed at 
right frequencies. 

We consider the experimental detection of these singu- 
larities, particularly fine structure effects, as the strong 
evidence in favor of the magnetic scenario for supercon- 
ductivity in the cuprates. 
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